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The joint extraction of iris and Vus from hadronic r decay data 
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We discuss the simultaneous determination of nis and Vus from flavor-breaking hadronic r decay sum rules, 
focussing on weight choices designed to better control problems associated with the slow convergence of the 
relevant integrated D — 2 OPE series. The results are found to display improved stability and consistency 
relative to those of conventional analyses based on the "(fc, 0) spectral weights". The results for rris agree well 
with those of recent strange scalar sum rule and strange pseudoscalar sum rule and lattice analyses. Results for 
Vus agree within errors with those of lattice-input-based r[Kfj,2]/T[-Kfj,2] and Kea analyses. Very significant error 
reductions are shown to be expected, especially for Vus, once improved strange spectral data from the B factory 
experiments becomes available. 



1. BACKGROUND 

Measurements of the inclusive hadronic r de- 
cay distributions for processes mediated by the 
flavor ij = ud, us vector (V) or axial vector 
(A) currents yield kinematically weighted linear 
combinations of the spectral functions, Py^j^.^j, 
of the spin J — and 1 parts of the rele- 



vant current-current correlators, 11 



{J) 



V/A-i] ■ 

t/^ hadronsy/A;y (7)]/r[T 
lyrE^Veij)], one has, explicitly [1], 



With 



Rv/A-ij = rfr 
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127r2|F,,|2W] 



dVr (1 - Vr) 



(1 + 2y.) p^Z^{s) 2yJ^]^.^^^{s)\ (1) 



where j/t — s/m^, Vij is the flavor ij CKM ma- 
trix element, Sew is a short-distance electroweak 
correction, and the superscript (0-1-1) denotes the 
sum of J = and J = 1 contributions. The ab- 
sence of kinematic singularities in the correlators 
corresponding to the spectral function combina- 
tions in Eq. (P) allows the spectral integrals to be 
re-written using the basic finite energy sum rule 



(FESR) relation. For such correlators H, with as- 
sociated spectral functions p, and w{s) analytic, 
this relation has the form 

/ dsw{s)p{s) = -1 / dsw{s)n{s) . (2) 
Jo 27rz J|s|=so 

Analogous FESR's, with spectral integral sides 
denoted Ryj^^^j, obtained from Rv/A;ij by rescal- 
ing the integrand with (1 — yr)'^y'!^ before in- 
tegration, are called the "(fc, m) spectral weight 
sum rules" . Similar spectral integrals and sum 
rules can be constructed for general non-spectral 
weights w{s), for either T^y^A^jis) or snj^j^. (s), 
and for any sq < m^. The corresponding spec- 
tral integrals are denoted generically by i?™ (so) 
in what follows. For "inclusive" sum rules (those 
with both J = and J = + 1 spectral contri- 
butions) the purely J = contribution will be 
referred to as "longitudinal" . 

Vus and/or rus are then to be extracted using 
flavor-breaking differences, SR'^{so), defined by 



SR'^iso) = 



(3) 
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Since SR"'{so) vanishes in the SU{3)p limit, its 
OPE representation, SRqp^{so), begins at di- 
mension D — 2. The D = 2 term is propor- 



tional to 



Experimental determinations of 



1 



2 



(5i?'"(so) over a range of Sq and w{s) allow a fit 
for TUs and/or V^s to be performed, so long as 
So is large enough that the OPE representation 
can be reliably employed |2l3l4l5l6l7j . The ap- 
proach is especially well-suited to the determina- 
tion of Vus :5 7^ since the smallness of makes 
the integrated D = 2 OPE contributions at scales 
'--^ 2 — 3 GeV^, and hence also the corresponding 
flavor-breaking spectral integral differences, much 
smaller than the individual flavor ud and us spec- 
tral integral terms (typically at the few to several 
percent level). Explicitly, since Eq. ([3|) implies 



IK., I = 



SR 



OPE 



(so) 



(4) 



it follows that an uncertainty A {SRqp^{so)) in 
SRqp^{sq) produces a fractional uncertainty ~ 
A {SRZpEi^o)) /2RZ{so) in IK. I which is much 
smaller than that on SRqp^{so) itself [5]. Moder- 
ate precision for 6Rqp^{sq) thus sufhces for high 
precision on |Ks|, provided experimental errors 
can be brought under control. 

We report here on a combined extraction of 
iris and |Ks| employing existing spectral data. 
The ud data, especially the V-l-A sum, are al- 
ready quite precise [8 9 lO llJ. The us data [TOl 
I12|13|14] . however, suffer from low statistics and 
have very sizeable errors above the K*. Ongoing 
analyses at BABAR and BELLE will very signif- 
icantly reduce the us errors in the near future. 
The focus here will be on bringing uncertainties 
on the theoretical (OPE) side of the analysis un- 
der better control, in particular those associated 
with the slower-than-previously-anticipated con- 
vergence of the relevant D = 2 OPE series [l5j . 

2. OPE COMPLICATIONS 

The first major problem on the OPE side is 
the very bad behavior of the integrated longti- 
tudinal D — 2 OPE series which shows no 
sign of converging at any kinematically accessi- 
ble scale [15116) , and, even worse, for all trunca- 
tion schemes used in the literature, badly violates 
spectral-positivity-based constraints [4ll7j . In- 
clusive analyses employing the longitudinal OPE 
representation are thus untenable. 



Fortunately, for a combination of chiral and 
kinematic reasons, the longitudinal D = 2 OPE 
problem is easily handled phenomenologically. 
Apart from the tt and K pole terms, longitudinal 
spectral contributions vanish in the SU{3)f limit 
and are doubly-chirally suppressed away from it. 
This double chiral suppression is preserved in 
the ratio of integrated non-pole to pole contribu- 
tions because the longitudinal kinematic weight 
in Eq. ^ has essentially the same value at the 
K pole as in the region of excited strange scalar 
and pseudoscalar (PS) resonances. The small 
residual non-pole us PS and scalar contributions 
can, moreover, be well-constrained phenomeno- 
logically [3], the former via a sum rule analy- 
sis of the us PS channel [18], the latter by dis- 
persive single- or coupled-channel -fCTr-scattering- 
data-based dispersive analyses [19|20j . With the 
dominant tt and K pole contributions already 
accurately known, a bin-by-bin subtraction of 
the longitudinal contributions to the experimen- 
tal distribution, and hence a determination of the 
(0 -|- 1) spectral function, can be performed, al- 
lowing FESR's not afflicted by the longitudinal 
D = 2 OPE problem to be constructed. 

Since the us scalar and PS spectral "models" 
used for the longitudinal subtraction correspond, 
via scalar and PS sum rules [18|20|21] . to values of 
TOs in excellent agreement with recent Nf = 2 -|- 1 
lattice results [22], significantly larger non-pole 
contributions are ruled out. The residual non- 
pole longitudinal subtractions are thus certainly 
small, and very well under control at the level 
required for (0 -I- 1) FESR analyses. We thus fo- 
cus, in what follows, on sum rules involving the 
flavor-breaking combination 



(o-Hi) 



(5) 



The D = 2 OPE contribution to AH is known 
up to 0[a^g) and given by [E] 



[An(Q^C"f = i9^[l + 2.333a-H 



19.9330^ + 208.7460^ + (2378 ± 200)a^] (6) 

where d = as((5^)/7r and = ms{Q^), with 
asiQ^) and ms{Q^) the running coupling and 
strange quark mass in the MS scheme. The 
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0(0**) coefficient is an estimate obtained using 
approaches which accurately predicted the 0{a^) 
coefficient in Eq. ^ and ny-dependent 0{a'^m^) 
coefficients of the electromagnetic current corre- 
lator in advance of their explicit calculation |23| . 
Since independent high-scale determinations of 
as{Mz) [21] correspond, after 4- loop running and 
matching |25|26j . to d{ml) ~ 0.10-0.11, the con- 
vergence of this series, at the spacelike point on 
the contour |s| — sq, is marginal at best, even at 
the highest scales accessible in hadronic t decay. 
Although |q;s((3^)| does decrease as one moves 
around the contour away from the spacelike point, 
allowing improvement in the convergence through 
judicious choices of weight, w(s), one must ex- 
pect to find very slow convergence of the inte- 
grated D = 2 series for those ^(s) not chosen 
specifically with this improvement criterion in 
mind. The {k,0) spectral weights, — 
(l-|-22/){l — 2/)'^+^, with y = s/sq, are highly non- 
optimal in this regard, since |1 — ?/| = 2\sin{(j)/2)\ 
(with 4> the angular position measured counter- 
clockwise from the timelike point) is peaked pre- 
cisely in the spacelike direction. Slow conver- 
gence, deteriorating with increasing k, is thus 
expected for the integrated D = 2 series of the 
(A;, 0) spectral weights. The results of Table I of 
Ref. [IS] bear out this expectation [27]. 

Several estimates of the integrated (0+1) D = 
2 OPE truncation uncertainty have been consid- 
ered in the literature: the size of the last term 
kept, the level of residual scale dependence, the 
difference between the truncated correlator and 
Adler function evaluations, and various combina- 
tions thereof. The slow convergence of the inte- 
grated series can make it difficult to be sufficiently 
conservative. E.g., the quadrature sum of the 
last term size plus residual scale dependence ver- 
sion of the O(a^) Adler function truncation uncer- 
tainty [5], yields a result ~ 2.5 times smaller than 
the actual difference between the 0(a'^)-truncated 
Adler function and 0(a*)-truncated correlator re- 
sults [7]. 

Since, due to the growth of ag with decreas- 
ing scale, higher order terms are relatively more 
important at lower scales, premature trunca- 
tion of a slowly converging series will induce 
an unphysical sg-dependence in extracted, nom- 



inally SQ-independent quantities. With polyno- 
mial weights, w(y) = l^m '"^ml/"^ (for which in- 
tegrated D — 2N + 2 OPE contributions not 
suppressed by additional powers of as scale as 
cn/sq) such unphysical SQ-dependence can also 
result if higher D contributions which might in 
principle be present are incorrectly assumed neg- 
ligible and omitted from the analysis. The ab- 
sence of phenomenological input for the values 
of the relevant D > 6 condensates makes such 
omission most dangerous for those w(y) having 
large coefficients c„i, with m > 2, where such 
unknown D > 6 contributions are potentially 
enhanced. The (2,0), (3,0) and (4,0) spectral 
weights, w(2:0)(-y) ^ I - 2y ~2y'^ + 8y^ - 7y^ + 2y'' , 
u,(3.o) (y) = 1 _ 3y + I0y3 _ i5y4 ^ _ 2y^ , and 

lly® + 2y'^ provide examples of weights having 
such large higher order coefficients. 

In view of the above discussion, SQ-stability 
tests are crucial to establishing the reliability of 
any theoretical error estimate. Failure to find a 
stability window in sq or, if not a stability win- 
dow, then at least a window within which the 
observed instability is safely smaller than the es- 
timated theoretical uncertainty, is a clear sign of 
an insufficiently conservative error. 

3. RESULTS AND DISCUSSION 

We restrict our attention, in what follows, to 
the V-l-A spectral combination, to weights satis- 
fying w(s — So) = 0, and to scales so > 2 GeV^, 
all of which serve to strongly suppress possible 
residual OPE breakdown effects 6,28,29 . 

The form of the known 13 = 4, 6 contributions 
to [An(Q^)]Qp^ may be found in Ref. [1]. Stan- 
dard values for the required OPE input param- 
eters are employed (for details, see Ref. [30]). 
Our central D = 2 determinations employ the 
contour-improved (CIPT) prescription [31] for 
the RG-improved correlator difference [An]jj^2 • 
An alternate CIPT evaluation, employing the 
truncated, RG-improved Adler function, provides 
one measure of the truncation uncertainty. (The 
two versions are equal to all orders but differ at 
0{a'^'^^) and higher when both are truncated at 
0{a^).) Exact solutions corresponding to the 4- 
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loop-truncated (3 and 7 functions [55] are used for 
the running of a and fhs- 

For the spectral integrals we employ the 
ALEPH ud [9] and us [T2] data, for which both 
data and covariance matrices are publicly avail- 
able. A small global renormalization of the ud 
data is required to reflect minor changes in the e, 
jjL and total strange branching fractions since the 
original ALEPH publication. Global mode-by- 
mode rescalings of the strange exclusive distribu- 
tions have also been performed to bring the orig- 
inal ALEPH branching fractions into agreement 
with current world average (PDG06) values [32] . 
Errors on the K and tt pole contributions have 
also been reduced by using the more precise val- 
ues implied by r[7r^2] and T[K^2]- The resulting 
errors on the ud and us spectral integrals are at 
the ^ 0.5% and 3 — 4% levels, respectively. 
The us errors will be drastically reduced by re- 
sults from BABAR and BELLE. 

3.1. The (fc,0) Spectral Weight Analyses 

As argued above, very slow convergence is ex- 
pected for the integrated D — 2, J = + 1 (fc, 0) 
spectral weight OPE series. Such slow conver- 
gence is seen explicitly in the results reported 
in Refs. [7115] . The situation for fc = is of 
particular practical interest since the sq = m^, 
(0, 0) us spectral integral is fixed by the total 
strangeness branching fraction. Spectral integral 
errors can thus be reduced through improvements 
to the various exclusive strange branching frac- 
tions. Such improvements are much less experi- 
mentally challenging than would be an improved 
determination of the full us spectral distribu- 
tion, which analyses employing other so and/or 
other weights would require. The very small nor- 
malization of the (0, 0) D — 2 OPE integral, 
which makes for a rather weak dependence of 
Vus on nis, is another favorable feature of this 
weight [5]. Unfortunately, these positive features 
must be weighed against the very poor conver- 
gence of the integrated D = 2 OPE series, and 
the concomitant difficulty in obtaining a reliable 
estimate of the truncation uncertainty. In the 
most recent version of this analysis (see the last 
of Refs. [S]) Vus is obtained using external input 
for iris, and sq = only. The quoted com- 



bined OPE-induced uncertainty is ±0.0011. As 
seen in Ref. [3, however, the (0, 0) sum rule dis- 
plays rather poor so-stability. Vus, e.g., changes 
by 0.0021 even over the rather restricted range 
— 0.4 GeV^ < sq < rri^. An alternate esti- 
mate of the D — 2 truncation uncertainty, ob- 
tained by taking twice the quadrature sum of the 
last term kept and the correlator-minus- Adler- 
fmiction difference yields a result, ±0.0022, much 
more in keeping with the size of the observed sq 
instability. Unfortunately, since a reduction in 
the truncation uncertainty does not appear likely, 
this more conservative estimate puts a sub-1% 
Vus determination out of reach of the (0, 0) spec- 
tral weight analysis. Figure 1 shows the OPE and 
spectral integral differences, as a function of sqi 
for various input = TOs(2 GeV), with Vus in 
each case obtained by fitting to the spectral in- 
tegral difference at sq = rn^. The very different 
SQ-dependences of the two curves is the source 
of the SQ-instability in the extracted Vus values. 
The figure makes clear that the instability in Vus 
found in Ref. [7] is not specific to the TOs employed 
as input in that analysis. 

A final illustration of the problematic features 
of the (fc, 0) spectral weight analyses is provided 
in Figure 4. Given any pair of (fc, 0) weights, rus 
and Vus can be obtained by a simultaneous fit us- 
ing the two So = spectral integrals as input. 
The figure shows the la contours for a series of 
such fits. It is evident that no good common fit 
region for rUs and Vus exists, further strengthen- 
ing the conclusion that the OPE representations 
for the (fc,0) spectral weights are not under suf- 
ficiently good control to allow a reliable determi- 
nation of nis and Vus ■ 

3.2. Non-Spectral Weight Analyses 

To reduce the problems associated with the 
slow convergence of the D = 2, J = 0±1 OPE 
series, we shift to FESR's based on three non- 
spectral weights, wio, ww, and W20, constructed 
in Ref. [3]. By design, these weights simultane- 
ously (i) improve D = 2 convergence; (ii) sup- 
press spectral contributions above 1 GeV^ (where 
us errors are large); and (iii) control weight co- 
efficients Cm, m > 2 (which might otherwise en- 
hance I? > 6 OPE contributions). Ref. [3J and 
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Table I of Ref. [7j show explicitly the much im- 
proved D = 2 convergence which results. Table II 
of Ref. [7] (the left, "AGO" , half corresponding to 
the us data treatment used here) also shows the 
much improved stability of Vus with respect to sq 
obtained for these weights, at least for the PDG04 
input TOs(2 GeV) = 105±25 MeV employed there. 
Figures 2 and 3 show the ms-dependent OPE vs. 
spectral integral comparisons, analogous to those 
shown for the (0,0) spectral weight in Figure 1, 
for W20 and wio, respectively. The results for wio, 
which are similar, have been omitted for brevity, 
but may be found in Ref. f30^. The existence 
of a window of values over which improved 
so-stability for Vus will be obtained for all three 
non-spectral weights is clear. Figure 5 shows the 
la contours for the various pairwise fits and com- 
bined 3-fold fit involving these weights. A good 
common fit region exists, in sharp contrast to 
the situation for the (fc,0) spectral weights. This 
strongly suggests that the lack of a good common 
fit region for the {k, 0) spectral weights is a result 
of the poor convergence behavior of the relevant 
integrated OPE series. 

The results of Figure 5 correspond to Vus = 
0.2202 ± 0.0046 and rris = ms(2 GeV) = 89 ± 25 
MeV. A modest reduction of the combined fit 
region can be achieved by adding a new non- 
spectral weight, wg,{y), to the analysis ^30J. An 
excellent common fit region remains. The cen- 
tral fit values and (somewhat) reduced errors are 
then [30] 

Vus = 0.2213 ± 0.0039, = 97 ± 19 MeV . (7) 

The result for Vus is compatible, within er- 
rors, with both (i) the IGIIEP06 Kis review up- 
date [33], = 0.2232 ±0.0006 (0.2249 ± 0.0019 
if the most recent lattice input for /+(0) [34] is re- 
placed by the old Leutwyler-Roos value [35]), and 
(ii) the recent MILG update of the r[is:^2]/r[7r^2] 
determination [36], Vus = 0.222ito"g°olt- The re- 
sult for TUs is in excellent agreement with recent 
strange scalar and PS sumrule, and Nf = 2-1-1 
lattice, results. 

To see how a reduction in the us spectral errors 
might impact the errors on Vus and TOs, we con- 
sider a scenario in which the us spectral function 
central values remain unchanged but the errors 



(covariances) are reduced by a factor of 3 (9). The 
combined fit errors on Vus and rus are reduced 
to ±0.0015 and ±13 MeV, respectively. This, of 
course, provides only a rough guide, since non- 
trivial shifts in the central us spectral function 
values are certainly to be expected. Nonetheless, 
the exercise makes clear that very significant re- 
ductions in the Vus error should be anticipated 
once B factory data becomes available. A more 
modest reduction is observed for the errors. 
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